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ABSTRACT 

We report the results of a survey of radio continuum emission of Galactic Wolf- 
Rayet stars north of 5 = —46°. The observations were obtained at 8.46 GHz (3.6 cm) 
using the Very Large Array ( VL A) , with an angular resolution of ~ 6" x 9" and typical 
rms noise of ~0.04mJybeam^^. Our survey of 34 WR stars resulted in 15 definite 
and 5 probable detections, 13 of these for the first time at radio wavelengths. All 
detections are unresolved {0^5"). Time variations in flux are confirmed in the cases 
of WR98a, WR104, WR105, and WR125. WR79a and WR89 are also variable in 
flux and we suspect they are also non-thermal emitters. Thus, of our sample 20 - 30 % 
of the detected stars are non-thermal emiters. Average mass loss rates determinations 
obtained excluding definite and suspected non-thermal cases give similar values for 
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WN (all subtypes) and WC5-7 stars (M(WN) = [4 ± 3] x lO'^ yr-^ and M(WC5- 
7) = [4 ± 2] X 10~^ Mq yr-i), while a lower value was obtained for WC8-9 stars (M(WC8- 
9) = [2 lb 1] X 10~^ Mq yr~^). Uncertainties in stellar distances largely contribute to the 
observed scatter in mass loss rates. Upper limits to the mass loss rates were obtained in 
cases of undetected sources or for sources which probably show additional non-thermal 
emission. 

Subject headings: radio continuum: stars; — stars: mass loss; — stars: Wolf-Rayet 

1. Introduction 

Wolf-Rayet (WR) stars are the evolved descendants of massive 0-type progenitors (Mi > 22 Mq) 
and are considered to be the immediate precursors of supernovae Type lb and Type Ic (e.g., Maeder 
1981; Meynet &i Maeder 2003). They arc characterized by powerful stellar winds with terminal 
velocities Uoo — 400-5500kms~-'^ (e.g., Prinja et al. 1990; Eenens & Williams 1994; Kingsburgh et 
al. 1995; for an overview see van der Hucht 2001, table 15, column 14), and mass loss rates in the 
range (1 - 5) x IQ-^ Mq yr^^ (e.g., Abbott et al. 1986 [hereafter AB86]; Leitherer et al. 1995, 1997 
[hereafter LC95, LC97]; Chapman et al. 1999 [hereafter CL99]). Due to these properties WR stars 
arc not only interesting in themselves, but also important species for chemical enrichment and 
galactic evolution (e.g., Dray et al. 2003; De Donder & Vanbeveren 2003). For general reviews of 
the WR phenomenon, see, e.g., van der Hucht (1992) and Maeder &; Conti (1994). 

The strong stellar winds of WR stars reveal themselves in the electromagnetic spectrum from 
the UV (P-Cygni profiles of resonance lines) to the IR/radio range (free-free radiation). Therefore, 
an important diagnostic tool to obtain information about stellar wind properties of massive stars is 
the study of the radio continuum emission from their close environs. Radio observations provide a 
reliable method for determining stellar mass loss rates if the emission is thermal in origin, since the 
measured radio flux density is related to the mass loss rate by an analytical expression (e.g., Wright 
&; Barlow 1975; Panagia &: Felli 1975). This emission is the free- free radiation originating in the 
ionized wind, expanding at constant velocity. It arises from the outer parts of the wind envelope 
and is partially optically thick in the radio region, showing positive spectral indices (a > -1-0.6, 
where 5,^ oc i/"). 

However, the last decade has seen an increasing number of discoveries of WR stars with non- 
thermal radio emission in addition to their free-free emission. Non-thermal emission is thought to 
be prodTiccd by electrons accelerated in strong shocks in the winds of single stars (White 1985) or in 
the shocked wind-wind collision region between massive WR+O binaries (e.g., Williams et al. 1990; 
Eichler & Usov 1993). Unlike thermal emission, non-thermal emission has a negative spectral index 
and can be variable. Van der Hucht et al. (1992) suggested that all cases of WR non-thermal radio 
emission correspond to WR-I-OB colliding wind binaries, which has been corroborated by Dougherty 
& WilUams (2000). 
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Due to their relatively large distances (van der Hucht 2001), the number of WR radio continuum 

detections was limited to only about 44 cases (AB86, LC95; LC97 and CL99). The VLA survey by 
AB86, which also quotes previous observations, contains radio data for 40 WR stars at 6> —47° 
with 27 detections using the VLA, with rms noise values of 0.06 - 0.23 mJy beam"^ . They found 77% 
probable or definite thermal wind sources and 12% probable or definite non-thermal sources. In the 
southern hemisphere, LC95, LC97 and CL99 surveyed 36 WR stars (S < —45°) at two frequencies 
using the ATCA, with rms noise values of 0.10-0.15mJybeam~^. Other WR radio continuum 
studies have been performed for WR 140 by Wilhams et al. (1987, 1990, 1994) and White & Becker 
(1995); for WR 146 by Dougherty et al. (1996, 2000) and Setia Gunawan et al. (2000); for WR147 
by Moran et al. (1989), ChurchweU et al. (1992), Davis et al. (1996), Wilhams et al. (1997), Skinner 
et al. (1999), Setia Gunawan et al. (2001), Watson et al. (2002) and Dougherty et al. (2003); for 
Galactic Center WR stars by Lang et al. (2001) and Lang (2003); for WR98a, WR 104 and WR 112 
by Monnier et al. (2002); and for massive stars in CygOB2 by Setia Gunawan et al. (2003b). In 
these studies, some 70 out of the ~250 known Galactic WR stars (van der Hucht 2001, 2003) have 
been investigated for radio continuum emission, resulting into 44 positive detections. The ATCA 
survey also indicates that at least 40% of the WR stars with measured radio spectral indices contain 
a non-thermal component at centimeter wavelengths, a percentage that, coincidentally, is equal to 
the percentage of WR-I-OB binaries among the known Galactic WR stars (van der Hucht 2001). 

In this paper we report on the results of a survey of radio continuum emission from Galactic 
WR stars at 3.6 cm with higher sensitivity than available before. Our aim is to increase the sample 
of WR stars with radio continuum observations in order to derive mass loss rates or upper limits, 
and to identify non-thermal cases. As a first step, we performed observations at one frequency only. 
Consequently, spectral indices of detected sources can only be determined here using radio data 
at other wavelengths from the literature and may be influenced by variability. Thus, new radio 
continuum observations are important to find new thermal and non-thermal WR radio sources, 
the thermal ones providing mass loss rates, crucial for testing existing wind models, and the non- 
thermal ones pointing to colliding winds in WR-I-OB binaries. 



2. Selected Wolf-Rayet stars 

In order to achieve these goals, we selected 33 WR stars from the catalog by van der Hucht 
(2001). To select the targets, we took into account all WR stars north of 5 =—46°, and selected the 
ones with expected free-free flux densities Scxp > 0.075 mJy. ^cxp was determined with the classical 
expressions by Panagia &: Felli (1975) and Wright & Barlow (1975) for the radio emission from 
an optically thick stellar wind flowing at constant velocity. In order to estimate these values, we 
adopt for all sources the same mass loss rate M = (1.5- 2.0)xl0~^ yr~^, and stellar distances 
and terminal velocities from van der Hucht (2001). 

With these assumptions, we found S^^p 

> 0.075 mJy for 64 WR stars. Only 27 of these stars 
had previous radio detections, and radio flux density upper Umits were available for 12 stars. We 
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were then left with 37 stars which either had not been observed yet or had not been detected. We 
selected 26 of these stars as our targets, listed in Table 1. In our sample we included some WR 
stars previously detected: WR79a (HD 152804, detected with the VLA by Bieging et al. 1989), 
WR81, WR89, WR98a, WR 104, WR105, and WR125. The detection of WR81 by AB86 was 
doubtful. WR98a, WR104, WR105 and WR125 are known non-thermal emitters and, thus, may 
vary in flux density. 



TABLE 1 



In order to find new non-thermal candidates, we included binaries and suspected binaries listed 
by van der Hucht et al. (2001). About 70% of the stars are closer than 3.0 kpc, while no stars beyond 
5.3 kpc are included. 

3. VLA observations 

We performed radio continuum observations at 8.46 GHz with the VLA in the direction of the 
selected WR stars. Thirty-three fields were observed on 15 September 2001 (during the move from 

the C- to the DnC-array), 5 and 8 October 2001 (in the DnC-array) and 12 November 2001 (in 
the D-array). With the exception of WR89, the field centers coincide with the stellar positions. 
The bandwidth was 50 MHz. The on-source observing time was about 25 minutes for each field to 
achieve a rms noise of ~0.04mJybeam~^. The synthesized beam is about 6'!lx&!0 in all cases. 

The data were edited, calibrated and imaged following a standard way using AiPS tasks. After 
editing, the data were calibrated both in amplitude and phase. The primary flux density calibrators 
were 3C 48 and 3C 286. We used phase referencing taking into account a selection of nearby 
point sources with precision positions (O^'l or better). The secondary calibrators were 0244-1-624, 
0747-331, 0828-375, 1001-446, 1604-446, 1744-312, 1751-253, 1820-254, 1832-105, 1851+005, 
1924+156, 2015+371, 2230+697 and 2250+558. The fields corresponding to WR4, WR96 and 
WR 142 were also self-calibrated in phase, since peak flux densities within the primary beam exceed 
lOmJybeam"^. Primary beam correction was applied to the images corresponding to WR89 and 
WR105. WR89 is in the same field as WR87, 1'.9 far from the field center. The primary beam 
correction factor for WR89 is 1.5. In the case of WR105, two additional point-like sources, with 
no relation to the WR star, are detected close to the center of the field. Diffuse emission from 
extended sources detected in the fields of WR87, WR104, WR106, WR113, WR114, WR121, 
WR 125, WR 142, WR 143 and WR 153ab was removed by applying an uv-range > 5KA (rejecting 
spatial scales ~40"). 
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4. Observed radio emission 

We report 15 definite and 5 probable detections at 3.6 cm out of the 34 observed WR stars. 
For 10 out of the 15 definite detections and 3 out of the 5 probable detections this represents a first 
radio detection, while four definite detections had been previously identified at 3.6 cm and one had 
been detected previously at a different frequency. All 20 detections arc unresolved. The images of 
the certain and probable detections are displayed in Figs. 1 to 20. The synthesized beam of each 
field is shown in a corner of the contour plots. 

Radio detections with signal-to-noise S / N > 5, as derived from the Gaussian fitting, and close 
agreement between radio and optical positions were classified as definite detections, while probable 
detections correspond to sources with S / N « 4 (WR 81, WR 95, WR 114 and WR 155) or with rela- 
tively large difference between the optical and the radio positions (WR79a and WR95). Detections 
classified as probable need further confirmation. 



FIGURES 1 - 20 



Flux densities and radio positions were obtained by fitting a 2-D Gaussian with the aips-jmfit 
tool. In most of the detections, the uncertainty in the radio position is less than 1". Except for 
WR79a and WR95, the optical and radio positions agree within 2". The flux density uncertainties 
quoted in the table are equal to the rms noise (Icr) of the images. Flux density uncertainties were 
derived from the images by taking into account a relatively large region near the detected sources. 
For the 14 undetected sources, the upper limits in flux density correspond to 3 a. Typically, the 
1 a level is in the range 0.025 - 0.084 mJybeam"^. In the cases of WR79a and WR89, the 1 a level 
is 0.13mJybeam~^. The large rms noise associated with the detection of WR89 is due to the fact 
that this star is displaced about 1^9 from the field center. The typical uncertainty in the radio 
position is 0.5". 

Table 2 lists the observed 3.6 cm flux densities and radio coordinates of the definite and probable 
detections, along with the optical coordinates from the catalog of van der Hucht (2001). Non- 
detected sources arc also indicated at the botton of the table. The latter have typical uncertainties 
of about 1" (Leitherer et al. 1997). 



TABLE 2 



In Table 3, we compare our 3.6 cm detections with previously determined flux densities or upper 
limits at six radio wavelengths. The table also gives spectral indices a when calculable. 
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Since we performed observations at one frequency only, spectral indices can not be derived 
from our data. Thus, the nature of the emission was inferred by comparing the observed fluxes 
with earlier determinations, allowing us to identify variable sources. Variability in flux density, 
which is a hint for the presence of non-thermal emission, is also indicated in Table 3. 



TABLE 3 



5. Comments on individual stars 

WR 79a, WN9ha (Fig. 4). This object is a known visual binary with Ac/) = 6'.'9 (Mason et al. 1998; 
van dcr Hucht 2001, table 22). The radio source extends to the W. The earliest flux density 
determinations by BA89 (VLA) indicated thermal emission with a spectral index a = +0.8. Recent 
ATCA observations (November 2000) by Setia Gunawan et al. (2003a) show a fiat non-thermal 
energy distribution with spectral index a = —0.0 ± 0.1. Our 3.6 cm observation, less than a year 
later, is consistent with the ATCA observation and, with respect to BA89 observation, with a 
variable non-thermal component. 

WR81, WC9 (Fig. 5). Our VLA image shows a source with a complex structure coincident with 
the optical position of the star. We identified the source to the N as the WR star because of 
its proximity (2'.'5) to the optical position. The radio source detected by AB86 and identified by 
these authors as the WR star (see their Fig. 1) is 4'.'6 south of the optical position and closer to 
the southern extension of the radio source detected in the new image. A spectral index cannot be 
derived. 

WR89, WN8h+0B (Fig. 7). This source was detected in the image corresponding to WR87, 
1^9 displaced from the field center. Previous flux density determinations by LC95 were consistent 
with a thermal wind. Judging from our observation, the 3.6 cm flux density dropped some 30 % in 
recent years. WR89 is associated with the open cluster Havlen-Moffat No. 1 (C1715— 387, Vazquez 
& Baume 2001). 

WR98a, WC8-9vd+? (Fig. 9). In the infrared this object shows a expanding dust spiral, rotating 
with a period P = 565ib50d (Monnier et al. 1999), which is considered to be the binary period. 
Monnier et al. (2002) observed this object with the VLA in 1999.7 and 2000.2 and found the source 
to be non-thermal. Our 3.6 cm flux density determination is ~ 20 % below their result. Apparently, 
the non-thermal source is variable. 

WR 104, WC9d+B0.5V (Fig. 12). In the infrared this object shows a expanding dust spiral, rotat- 
ing with a period P = 243 it 3 d (Tuthill et al. 1999), which is considered to be the binary period. 
Monnier et al. (2002) found the source to be non-thermal. Our 3.6 cm flux density determination 
is ~ 40 % below theirs, thus, the non-thermal source is apparently variable. 
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WR105, WN9h (Fig. 13). This star shows variable non-thermal flux densities. Our 3.6 cm flux 
density is ~ 50 % larger than the one measured by LC97. The double radio source to the S of the 
optical position of WR 105 was also detected by AB86 (see their fig. 1). These radio sources were not 
detected in the NVSS (Condon et al. 1998), nor were optical counterparts found. Zoonematkermani 
et al. (1990, radio continuum at 20 cm) detected the double source with flux densities of 22 mJy and 
16 mJy, while from our VLA data we estimate 6.4 it 0.2 mJy and 5.4 it 0.2 mJy. These non-thermal 
sources have spectral indices a of —0.7 and —0.6, respectively. 

WR125, WC7ed+09III (Fig. 17). The AB86 1982 observations at 6.3 cm and 2.0 cm showed a 
non-thermal radio source with a = — 0.5. Our 2001 3.6 cm flux density observation agrees with the 
AB86 data for a non-thermal source with that spectral index. The source is variable in flux and 
turned from non-thermal to thermal in WH92 (see Table 3). We suggest here that the emission is 
again non-thermal because of the similarity in flux density levels between the AB86 and the present 
data. The detection of non-thermal emission in 1982 and in 2001 suggests a period of around 20 
years. On the other hand, the last IR excess of the WC-l-0 binary WR 125 started in 1991 (WH92), 
continued till 1993 and then faded away. IR monitoring performed during the last 18 years showed 
only one period of infrared excess (Williams 2002), implying an orbital period of at least 18 years. 
WR 125 has much in common with the archetype WR colliding wind binary WR 140 (HD 193793, 
WC7pd-|-04-5, Williams et al. 1990), which has a thermal to non-thermal radio light curve with 
a period of 7.94 yr (Williams et al. 1994). If WR125 is an analog of WR140, then its 1992 IR 
excess occurred at periastron passage in an eccentric orbit, and an IR excess will happen again in 
~2010-2012. 

WR133, WN5-I-09I (Fig. 18). We detect a double point-like source (separation ~ 10") with the 
southern component coincident with the optical position of the WN5-I-09I binary (Fig. 18, upper 

panel), which is a known visual binary (A(^ = 5'.'4, Hartkopf et al. 1999; van der Hucht 2001, table 
22). The source to the NE has a flux density of 6*3.60111 = 0.39 ± 0.03 mJy bcam^-*^. The lower panel of 
Fig. 18 displays an overlay of the radio continuum emission and the DSS2 red image, suggesting the 
presence of an additional star NE of the WR system, coincident in position with the northern radio 
source. Pollock et al. (1995) found a ROSAT X-ray excess for WR133, indicative for a colliding 
wind binary. 

Finally, we note that a radio source with a flux density S^ Qcm — 

0.15 ± 0.04 mJy is detected 2'.'5 

far from the optical position of WR4. We believe that this star has probably been barely detected, 
but, because of the low S / N ratio and the difference between the optical and radio positions (2.1cr), 
we have listed this star as non-detected. 

In summary, we note that WR79a, WR89, WR98a, WR104, WR105 and WR125 are vari- 
able, which can be explained by the presence of a non-thermal component. This implies that our 
observations of these six stars cannot really be used for mass loss rate determinations, or will yield 
only upper limits. Of the above six cases all three WCd objects (WR98a, WR104 and WR125) 
turn out to be variable non-thermal radio sources, suggesting a large number of non-thermal cases 
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between these objects. Additional observations are necessary to investigate the nature of the other 
three WCd stars in our sample (WR95, WR103 and WR113) since previous observations only 
provided upper limits to the flux densities. 



6. Radio mass loss rates 

Following Panagia & Felli (1975) and Wright & Barlow (1975), the flux density due to thermal 
Bremsstrahlung in an ionized expanding stellar wind can be expressed as 

= 2.32 X 10^(^^)4/3(1^)2/3 (1)^ 

where S,^ is the flux density in mJy, Af is the mass loss rates in MQyr~^, is the terminal 
velocity in kms^^, v is the frequency in Hz and d is the distance in kpc. fi, Z and 7 are the mean 
molecular weight, the rms ionic charge and the mean number of electrons per ion, respectively. This 
expression corresponds to a spherically symmetric, stationary, isothermal wind flowing at constant 
velocity. The free-free Gaunt factor gg can be approximated, following Leitherer &; Robert (1991) 
as 

rp3/2 

gs = 9.77(1 + O.lSlog^) (2). 

In this expression, Tg is the electron temperature of the wind in K. 

Based on Equations (1) and (2), we have derived mass loss rates M corresponding to the 
observed flux densities at 3.6 cm, assuming that the observed radio emission is due to free-free 
emission. In the case of non-thermal contributions, our values represent upper limits to the true 
mass loss rate. All non-detections are indicated as upper limits to M. Following LC97, we adopted 
Tg = lO^K in expression (2). As pointed out by these authors, this temperature has a very minor 
effect on the gg. Values for the mean molecular weight, the rms ionic charge and the mean number 
of electrons per ion were taken from LC97. The values adopted for the WO star in our sample, 
WR142, are those for the WC4 subtype (Crowther et al. 1998). 

Table 4 lists the derived mass loss rates M and upper limits. Uncertainties in these values come 
from several sources. Errors in mass loss rate determinations were estimated following eq. (5) in 
LC97. We adopted similar criteria as LC97 to estimate the uncertainty in each parameter. For the 
error in the distances we assumed 20% in case of cluster/associations distances and 50% in case 
of field stars. The logarithmic uncertainty is then ±0.09 and ±0.25, respectively. Errors in gg are 
about 10% (0.05 in logarithmic scale). We adopted uncertainties in wind terminal velocities of 10% 
(i.e., a logarithmic error of ±0.04). Finally, logarithmic errors in Z, 7 and n were adopted as ±0.08. 
Typical logarithmic errors in logM are in the range ± 0.20 to ± 0.24 in case of cluster /association 



-9- 



distances and ± 0.40 in case of field stars. With the quoted uncertainties for the stellar distances 
and other wind parameters, uncertainties in mass-loss rates determinations are dominated by errors 
in distances, while uncertainties in flux densities have a minor impact. 



TABLE 4 



Other sources of error are non-thermal contributions to the flux at 3.6 cm. Inhomogeneities 

or clumping in the stellar winds lead to over estimates of the final values as emphasized by, e.g., 
Contrcras ct al. (1996) and Morris ct al. (2000). The latter found volume filling factors of the range 
0.04 - 0.25, leading to mass-loss rates a factor of 2.5 - 5 lower than expected for smooth, homogeneous 
winds. These effects have not been considered here. 



7. Discussion 

In order to compare mass loss rates derived for different WR subtypes from available radio 
surveys, we plotted the ATCA and VLA mass loss rates versus WR subtype in Fig. 21. The ATCA 
surveys have been performed by LC95, LC97 and CL99. Previous VLA results were published by 
Willis (1991) (who updated AB86 results by taking into account new criteria in terminal velocities 
and in the chemical composition of the region where radio emission originates) and by Monnier et 
al. (2002). 



FIGURE 21 



Monnier ct al. observed WR 98a and WR 104 at several frequencies and found a non-thermal 
component in addition to the thermal one. They separated both components and derived a mass- 
loss rate from the best-fitting value to the thermal emission. The M-values they derived are 
0.8x10"^ Moyr-i for WR98a and 0.5x10"^ Moyr'^ for WR104. 

The upper panel of Fig. 21 displays the mass loss rate (log M) as a function of the spectral WR 
subtype for the stars included in our VLA survey. Bona-fide mass loss rate estimates are indicated 
by filled circles, while open circles indicate upper limits both for the variable detected sources and 
the undetected sources. Uncertainties in the derived values are also indicated. Note that we have 
plotted the mass loss rates derived by Monnier et al. instead of our estimates listed in Table 4. 

The lower panel of Fig. 21 shows all radio mass loss rates from the ATCA and VLA samples. 
Although more than one mass loss rate determination is available for some of the stars, we plotted 
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only one value for each star from (in order of preference) Monnier et al. (2002), new VLA data, 
ATCA data and previous VLA surveys. In the cases of non-thermal sources listed by LC97 and 
CL99, the M-valucs plotted arc upper limits. As for the upper panel, open circles indicate upper 
limits both for undetected sources and variable detected sources. 

Excluding definite and suspected non-thermal cases (WR 79a, WR 89 and WR 105) the average 
mass loss rate for the WN stars detected in our survey is M(WN) = [4.3 ± 2.5] x 10^^ Mq yr~^, equal 
to the average WN mass loss rate in ATCA and previous VLA data sets of M(WN) = 4x 10~^ M0 yr~^ 

Mass loss rates determined from emission lines studies by Nugis &: Lamers (2000) differ from 
values derived from radio continuum data: in some cases (WR104 and WR113), M derived from 
emission lines are factors of two smaller than our radio estimates, while in other cases (WR81, 
WR95 and WR103) factors of two larger than our results. Clumping and/or an ionization or 
temperature structure different from our assumptions may account for the differences. 

We have thirteen WC8-9 stars included in our sample and detected seven of them. Four of 

these stars had no radio detections yet (WR88, WR95, WR103 and WR113), thus increasing to 
12 the number of WC8-9 stars with radio detections. Most of these stars have heated circumstellar 
dust shells (sec Williams 1995). WR95, WR103 and WR113 arc known to have persistent dust 
shells (see Table 1). The distribution of these stars in the diagram displays the same trend as in 
fig. 17 of LC97, with relatively low mass loss rates. The average mass loss rate for WC8-9 stars 
detected in our sample, excluding definite non-thermal cases (WR98a and WR104) is M(WC8- 
9) = [1.7 ± 1.0] X 10~^ M0 yr~^ . The average mass loss rate for WC5-7 stars, obtained from bona-fide 
values from the ATCA (5 stars) and the previous (2 stars) and present (1 star) VLA samples, is 
M(WC5-7) = [4.4 ± 1.6] X 10-5 Moyr'^ higher than the value derived for WC8-9 stars. Note, on 
the other hand, that average values for WN and WC5-7 stars are equal. Uncertainties in stellar 
distances clearly contribute to the observed scatter of mass loss rates. 

The presence of dust in the wind envelopes might decrease the escaping UV photon flux, 
lowering the thermal radio emission, as suggested by Monnier et al. (2002), and thus the derived 
mass loss rates. LC97 suggest that this may explain the low mass loss rates obtained for WC9 stars 
from radio data. Note that the mass loss rates for the WC8 star WR98a, which also displays a 
dust persistent shell, is also relatively low (see Fig. 21, lower panel). 

The percentage of non-thermal emitters between the detected sources in our sample is 20- 
30%. A similar value is derived for the 1986 VLA survey (which includes also previous results). 
Our percentage is lower than the one derived from ATCA data (40%). 

The case of WR89 (suggested to be a binary by van der Hucht 2001), which shows time 
variations in flux density (see Table 3), hints to additional non-thermal emission (the radio spectrum 
was thermal at the time of the ATCA observations) . WR 79a also appears to be variable in flux, 
suggesting the presence of a non-thermal component. However, a variable thermal component can 
not be discarded (e.g., Watson et al. 2002). It would be desirable to perform multifrequency radio 
observations of WR 79a and WR 89 to look for variability and to analyze the nature of the flux 
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variations. 

Finally, we would like to comment on the undetected stars. We repeat here that uncertainties 
in the expected flux densities S'exp are quite large, arising mainly from uncertainties in distances 
of WR field stars and in the expected range of mass loss rates (1 - 5x 10~^ Mq yr~^). Scxp values 
for WR79b, WR92, WR96 and WR119 are 0.12-0.26, 0.09-0.14, 0.10-0.16 and 0.11 - 0.16 mJy, 
respectively, which are below our detection limit (see Table 2). ^exp values were estimated adopting 
the stellar parameters of Tables 1 and 4, and mass loss rates in the range (1.5 - 2.0) x 10~^ Mq yr~^. 
For the cases of WR80 and WR121, however, Sexp = 0.5 - 0.8 and 0.4-0.6 mJy, respectively. The 
derived upper limit is 0.17 mJy (Scr) for both stars. Although a mass loss rate lower than typical 
can not be discarded, the objects may be at larger distances than adopted here, thus lowering the 
expected flux density. 

8. Conclusions and suggestions for future work 

We performed a survey of 34 Galactic Wolf-Rayet stars at 3.6 cm within the declination range 
of the Very Large Array. We report 15 definite and 5 probable detections. Of these 20 sources, 13 
were detected for the first time at radio frequencies. 

We confirm time variations in fiux of WR98a, WR104, WR105, and WR125, which support 
the presence of a non-thermal component. For WR 125 (WC7ed-|-09III) we suggest a binary period 
lower than 20 yr. WR79a (HD 152408, WN9ha) and WR89 are also variable in flux and we suspect 
they are also non-thermal emitters. Non-thermal radiation is indicative for wind-wind collision 
in WR+OB binaries. Thus, of our sample 20-30% of the detected stars are non-thermal radio 
sources. 

Averages of mass loss rate determinations yield similar values for WN (all subtypes) and WC5- 
7 stars (M (WN) = [4 ± 3] x 10"^ Mq yr'^ and M(WC5-7) = [4 ± 2] x 10"^ Mq yr'^), while a slightly 
lower value was derived for WC8-9 stars (M(WC8-9) = [2 ± 1] x 10~^ M© yr"^). 

Future observations at several radio frequencies should monitor the objects WR79a, WR89, 
WR98a, WR104 and WR125 during long time periods. These observations can help to determine 
the presence of non-thermal contributions to the flux density as a function of time. Radio obser- 
vations at low and high frequencies of the sources first detected will be useful to investigate the 
nature of the radio emission. 
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Table 1: Selected Wolf-Rayet stars 



WR 


other 


spectral 


d 


binary 


P 


Voo 


rf 


(l(*si,;^iialj()ii (s) 


lyix- 


(kpc) 


si al us 


(d) 


(kms"'^) 


4 


HD 16523, V493 Per 


WC5+? 


2.44 


SBl, no d.e.l. 


2.4 


1900 


5 


HD 17638 


WC6 


1.91 






2100 


8 


HD 62910 


WN7/WCE+? 


3.47 


SBl 


38.4 


1590 


12 


Vc5-5 


WN8h+? 


5.04 


SBl, no d.c.h 


23.9 


1100 


79a 


HD 152408 


WN9ha 


1.99 


VB (6'.' 9) 




935 


79b 


HD 152386 


WN9ha 


2.90 


VB (0'.'55) 




1650 


80 


Wral581 


WC9d 


1.50 






1200 


81 


He3-1316 


WC9 


1.57 






910 


82 


LSll 


WN7(h) 


5.25 




2.0 


1100 


87 


LSS4064 


WN7(h)+0B 


2.88 


a, d.e.l. 




1400 


88 


Thcl 


WC9 


2.33 






1125 


89 


AS 223 


WN8h+OB 


2.88 


a, d.e.l., VB (9'.'88) 




1600 


92 


HD 157451 


WC9 


3.84 






1100 


95 


He3-1434 


WC9d 


2.09 






1100 


96 


LSS4265 


WC9d 


3.58 






1100 


98a 


IRAS 17380-3031 


WC8-9vd+? 


1.9 


CWB 


565 


900i 


100 


HDE 318139 


WN7 


4.61 






1600 


103 


HD 164270, V4072 Sgr 


WC9d+? 


2.21 


SBl, no d.e.l. 


1.8 


1100 


104 


Ve2-45 


WC9d+B0.5V 


2.3 


SB2, VB (0'.'975) 


243.5 


1220 


105 


Ve2-47 


WN9h 


1.58 






700 


106 


HDE 313643 


WC9d 


2.32 






1100 


113 


HD 168206, CVSer 


WC8d+08-9IV 


1.79 


SB2 


29.7 


1700 


114 


HD 169010 


WC5+OB? 


2.00 


d.e.l. 




2000 


119 


The 2 


WC9d 


3.31 






1200 


120 


Vyl-3 


WN7 


3.56 






1225 


121 


AS 320 


WC9d 


1.83 






1100 


124 


209 BAG, QRSge 


WN8h 


3.36 


SBl?, no d.e.l. 


2.7 


710 


125 


IC14-36, V378Vul 


WC7ed+09ni 


3.06 


SB2 


>6600 


2900 


133 


HD 190918, V1676 Cyg 


WN5+09I 


2.14 


SB2, VB (5^'4) 


112.4 


1800 


142 


Sand 5 


W02 


0.95 






5500^ 


143 


HD 195177 


WC4+OB? 


1.07 


d.e.l. 




2750 


153ab 


HD 211853, GPCcp 


WN6/WCE +06I 


2.75 




6.7 


1785 


155 


HD 214419, CQCep 


WN6+09II-Ib 


2.75 


SB2 


1.6 


1400 


156 


AC +60 38562 


WN8h+OB? 


3.56 


d.e.l. 


6.5 


660 



Notes: 

Cols.l and 2: WR catalogue number and any other designation(s) ; 
col. 3: spectral classification; col. 4: heliocentric distance; 

col. 5: the binary status, in case of visual binaries (VB), the angular separation is indicated; 
col. 6: period; 

col. 7: terminal wind velocity. 

a: absorption lines present in optical spectrum; 

d.e.l.: diluted emission lines 

All parameters and values from van der Hucht (2001) except: 
1: from Monnier et al. (2002); 
2: from Kingsburgh et al. (1995). 
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Table 2: Optical positions, radio positions, and flux densities of the observed Wolf-Rayet stars 



WR optical position radio position /Ss.ecm epoch 



# R.A.(J2000) Dec.(J2000) R.A.(J2000) Dec.(J2000) 







h m 


s 


O 




// 




h m s 




(mJy) 
























Definite detections 








5 


02 


52 


11, 


,66 


+56 


56 


07, 


.1 


02 52 11.62 + 0.05 


+56 56 08.2 + 0.4 


0.20 + 0.03 


2001, 


.9 


s 

O 


07 


AA 


PIS 

GO. 


22 




54 


29 




07 44 58.26 + 0.02 


-31 54 29.7 + 0.2 


n 4C _|_ n no 


2001 

ZiUU -L - 


s 

.O 


12 


08 


44 


47, 


,25 


-45 


58 


55, 


.5 


08 44 47.33 + 0.03 


-45 58 52.5 + 0.7 


0.51 + 0.06 


2001, 


.8 


88 


17 


18 


49, 


.50 


-33 


57 


39, 


.8 


17 18 49.67 + 0.05 


-33 57 40.9 + 0.6 


0.26 + 0.05 


2001, 


.8 


89 


17 


19 


00, 


.52 


—38 


A O 

48 


51, 


o 

.2 


17 19 00.49 + 0.03 


-38 48 49.4 + 0.4 


o r\ 1 r\ 1 

2.U +U.1 


2001, 


.7 


aoa 


1 '7 


41 


12, 


.9 


QO 

— oU 


QO 


on 
29 




17 41 13.08 + 0.03 


-30 32 30.0 + 0.3 


o /i'7_i_ o riK 
U.4/ ± U.Uo 


onm 
2UUi, 


Q 


iUU 


1 T 

17 


42 


09, 


,77 


—32 


o o 
33 


O A 

24 


.7 


17 42 09.78 + 0.02 


-32 33 24.9 + 0.2 


U.47 + U.04 


OA A 1 


o 
.0 


1 no 




01 


43, 


.14 


OO 

—6Z 


A O 


55, 


o 
.2 


18 01 43.20 + 0.05 


-32 42 55.7 + 0.7 


A 1 _l_ A A /t 

0.21 + 0.04 
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. 1 


104 


18 


02 


04, 


.07 


-23 


37 


41, 


.2 


18 02 04.16 + 0.02 


-23 37 41.0 + 0.3 


0.54 + 0.06 


2001, 


.8 


105 


18 


02 


23, 


,46 


-23 


34 


37, 


.7 


18 02 23.48 + 0.01 


-23 34 37.3 + 0.3 


5.4 + 0.1 


2001, 


.9 


113 


18 


19 


07, 


,36 


-11 


37 


59, 


.2 


18 19 07.40 + 0.01 


-11 37 58.9 + 0.1 


0.75 + 0.04 


2001, 


.8 


120 


18 


41 


00, 


.88 


-04 


26 


14, 


.3 


18 41 00.85 + 0.03 


-04 26 14.6 + 0.4 


0.40 + 0.04 


2001, 


.9 


125 


19 


28 


15, 


.57 


+19 


33 


21, 


.1 


19 28 15.61 + 0.01 


+19 33 21.6 + 0.1 


1.14 + 0.03 


2001, 


.9 


133 


20 


05 


57, 


,33 


+35 


47 


18, 


.2 


20 05 57.31 + 0.02 


+35 47 18.6 + 0.3 


0.36 + 0.03 


2001, 


.9 
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23 


00 


10, 


.13 


+60 


55 


38, 


.4 


23 00 10.14 + 0.01 


+60 55 38.4 + 0.1 


1.06 + 0.03 


2001, 


.9 




















Probable detections 








79a 


16 


54 


58, 


,51 


-41 


09 


03 


.1 


16 54 58.53 + 0.05 


-41 09 07.8 + 0.8 


0.68 + 0.13 


2001, 


.7 


81 


17 


02 


40, 


,39 


-45 


59 


15, 


.5 


17 02 40.48 + 0.06 


-45 59 12.9+ 1.4 


0.29 + 0.07 


2001, 


.8 


95 


17 


36 


19, 


.76 


-33 


26 


10, 


.9 


17 36 19.87 + 0.05 


-33 26 13.4 + 0.7 


0.16 + 0.04 


2001, 


.7 


114 


18 


23 


16, 


.39 


-13 


43 


25, 


.8 


18 23 16.42 + 0.04 


-13 43 25.9 + 0.2 


0.15 + 0.04 


2001, 


.8 


155 


22 


36 


53, 


.96 


+56 


54 


21, 


.0 


22 36 53.79 + 0.10 


+56 54 21.5 + 0.7 


0.12 + 0.03 


2001, 


.9 




















Undetected sources 








4 


02 


41 


11, 


.68 


+56 


43 


49, 


.7 






<0.15 


2001, 


.9 


79b 


16 


55 


06, 


.45 


-44 


59 


21, 


.4 






<0.20 


2001, 


.7 


80 


16 


59 


02, 


,2 


-45 


43 


06 








<0.17 


2001, 


.8 


82 


17 


04 


04, 


,61 


-45 


12 


15, 


.0 






<0.23 


2001 


.7 


87 


17 


18 


52, 


,89 


-38 


50 


04, 


.5 






<0.24 


2001, 


.7 


92 


17 


25 


23, 


.15 


-43 


29 


31, 


.9 






<0.18 


2001, 


.8 


96 


17 


36 


24, 


,2 


-32 


54 


29 








<0.16 


2001, 


.7 


106 


18 


04 


43, 


.66 


-21 


09 


30, 


.7 






<0.17 


2001, 


.8 


119 


18 


39 


17, 


.91 


-10 


05 


31, 


.1 






<0.11 


2001, 


.7 


121 


18 


44 


13, 


.15 


-03 


47 


57 


.8 






<0.17 


2001, 


.9 


124 


19 


11 


30, 


.88 


+16 


51 


38, 


.2 






<0.25 


2001, 


.9 


142 


20 


21 


44, 


.36 


+37 


22 


30, 


.3 






<0.90 


2001, 


.9 


143 


20 


28 


22, 


.68 


+38 


37 


18, 


.9 






<0.12 


2001, 


.9 


153ab 


22 


18 


45, 


.61 


+56 


07 


33, 


.9 






<0.14 


2001, 


.9 
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Table 3: Radio flux densities of earlier detected Wolf-Rayet stars, compared with our detections 



WR 




present and previous observations (mjy) 




spectral 




# 




5l2.6cin 


Se.Scm 


S3. 6cm 


S2.0cin 


■Si. 3cm 


index a 




79a 




1.0±0.1 


1.1±0.1 
0.8±0.1 


0.9 ±0.1 

0.7 ±0.1 


2.4±0.1 




+0.8 
-0.0 
V 


BA89 
SC03 
this study 


81 






0.3 ±0.08 


0.3 ±0.1 








AB86 
this study 


88 


<1.29 


<0.57 


<0.42 


<0.42 
0.26 ±0.05 








LC97 
CL99 
this study 


89 


<1.20 


<0.90 


0.6 ±0.1 
1.9 ±0.2 


3.0 ±0.1 
2.0 ±0.1 






+0.8 
V 


AB86 
LC95 
CL99 
this study 


95 


<1.50 


<0.66 


<0.4 
<0.45 


<0.45 
0.16 ±0.04 








AB86 
LC97 
CL99 
this study 


98a 


<0.36 




0.37 ±0.07 


0.60 ±0.05 
0.47 ±0.05 


0.64±0.11 


0.57 ±0.10 


+0.3 
V 


MT02 
this study 


103 


<0.90 


<0.42 


<0.2 
<0.42 


<0.42 
0.21 ±0.04 








AB86 
LC97 
CL99 
this study 


104 


<1.59 
<0.30 


<0.99 


<0.4 
<2.01 


<0.39 

0.87 ±0.06 
0.54 ±0.06 


1.02 ±0.12 


0.94 ±0.10 


+0.1 
V 


AB86 
LC97 
CL99 
MT02 
this study 


105 


<1.17 


<0.69 


3.6 ±0.21 
4.39 ±0.15 


3.8 ±0.2 
5.4 ±0.1 






-0.3 
V 


AB86 
LC97 
CL99 
this study 


113 


<2.25 


<0.90 


<0.4 
<0.80 


<0.80 
0.75 ±0.04 








BA82 
LC97 
CL99 
this study 


114 


<1.17 


<0.54 


<0.3 
<0.45 


<0.45 
0.15 ±0.03 








AB86 
LC97 
CL99 
this study 


125 


1.53 ±0.06 




1.5 ±0.09 
1.18 ±0.06 
0.20 ±0.04 


1.14 ±0.03 


0.8 ±0.09 
0.82 ±0.10 




-0.5 
-0.3 
+0.7 V 


AB86 
WH92 
WH92 
this study 


133 






<0.3 


().:«) ±0.0:-! 








AB86 

1 liis study 



Note: V: apparently variable. References: AB86: Abbott et al. (1986); BA82: Bieging et al. (1982); BA89: 
Bieging et al. (1989); CL99: Chapman et al. (1999); LC95: Leitherer et al. (1995); LC97: Leitherer et 
al. (1997); MT02: Monnier et al. (2002); SC03: Setia Gunawan et al. (2003a); WH92: Williams et al. (1992). 



Table 4: Mass loss rates from 3.6 cm observations of our sample of 34 Wolf-Rayet stars 
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d fx Z 
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logM 
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# 
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Fig. 1. — Contour image for WR5. The eross marks the optical position of the WR star, and the 
elUpse in one of the corners, the synthesized beam. We note that the extend of the cross does not 
represent the uncertainty in the stellar position. Contour lines are -3 (dashed contour), 3, 4, 5, 6 
and 7(7 (la = 0.025 mJybeam"^). 

Fig. 2. — Contour image for WR8. Contour lines are -3 (dashed contour), 3, 4, 5, 6, 8, 10, 12 and 
14cr (Icr = 0.030 mJybeam-^). 

Fig. 3. — Contour image for WR12. Contour lines are -3 (dashed contour), 3, 4, 5, 6, 7, 8 and 9a 
(1(7 = 0.055 mJybeam"^). 

Fig. 4. — Contour image for WR79a. Contour lines are -3 (dashed contour), 3 and 4(7 (Icr = 0.125 
mJy beam~^). 

Fig. 5. — Contour image for WR81. Contour lines are -3 (dashed contour) and 3c7 (Icr = 0.065 
mJy beam~^). 

Fig. 6. — Contour image for WR88. Contour lines are -3 (dashed contour), 3, 4 and 5cr (Icr = 
0.050 mJybeam"^). 

Fig. 7. — Contour image for WR89. Contour lines are -3 (dashed contour), 3, 4, 5, 6, 7, 9, 11 and 
13cr (Icr = 0.130 mJybeam-i). 

Fig. 8. — Contour image for WR95. Contour lines are -3 (dashed contour) and 3cr (Icr = 0.040 
mJy beam~^). 

Fig. 9. — Contour image for WR98a. Contour lines are -3 (dashed contour), 3, 4, 5, 6 and 7cr (Icr 
= 0.050 mJybeam"^). 

Fig. 10. — Contour image for WRIOO. Contour lines are -3 (dashed contour), 3, 4, 5, 6, 7, 8, 10 
and 12(7 {la = 0.040 mJybeam'^). 

Fig. 11. — Contour image for WR103. Contour lines are -3 (dashed contour), 3 and 4cr (Icr = 
0.040 mJybeam"^). 
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Fig. 12. — Contour image for WR104. Contour lines are -3 (dashed contour), 3, 4, 5, 6 and 7a 
{la = 0.060 mJy beam~^). 

Fig. 13. — Contour image for WR 105. Contour lines are -3 (dashed contour), 3, 6, 12, 24, 36, 48, 
and 60c7 (la = 0.080 mJybeam^i). 

Fig. 14. — Contour image for WR113. Contour lines are -3 (dashed contour), 3, 6, 9, 12, 15 and 
17(j (la = 0.040 mJybeam-i). 

Fig. 15. — Contour image for WR114. Contour lines are -3 (dashed contour) and 2>a {la = 0.040 
mJy beam~^). 

Fig. 16. — Contour image for WR120. Contour lines are -3 (dashed contour), 3, 4, 5, 6, 7, 8, 9 
and lOo" (la = 0.035 mJybeam^^). 

Fig. 17. — Contour image for WR125. Contour lines are -3 (dashed contour), 3, 6, 9, 12, 15, 18, 
24, 30 and 35a (la = 0.030 mJybeam^^). 

Fig. 18. — Contour image for WR133. Upper panel: Contour lines are -3 (dashed contour), 3, 4, 
6, 8, 10, 12 and 14a (la = 0.025 mJybeam"^). Lower panel: Overlay of the DSS2 red image of 
WR 133 (gray scale, in arbitrary units) and the radio continuum emission at 3.6 cm. 

Fig. 19. — Contour image for WR155. Contour lines are -3 (dashed contour), 3 and 4a (la = 
0.025 mJybeam-i). 

Fig. 20. — Contour image for WR156. Contour lines are -3 (dashed contour), 3, 6, 9, 12, 15, 18, 
24 and 30a (la = 0.030 mJybeam"^). 

Fig. 21. — WR mass loss rates (log M, where M is in units of MQyr~^) from radio observations 
versus spectral subtype. Upper panel: Mass loss rates derived from the present VLA survey. Bona 
fide mass loss rates values corresponding to detected sources are indicated as filled circles. Bar 
errors are included. Open circles indicate upper limits corresponding to both detected sources with 
variable flux density and undetected sources. Lower panel: Mass loss rates derived for the VLA 
and ATCA target stars. Open and filled circles have the same meaning as in the upper panel. 
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